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PREFACE 


The present thesis entitled "A study of some 
problems-j^of blood flow in human circulatiory system", 
is a record of genuine research work carried out by- 
me under the supervision and guidance of Dr. J.S. 
Parihar (Supervisor), Reader in the Deptt. of Maths, 
Atarra P.G. College Atarra (Banda) 

The purpose of the thesis is to study some 
problems on blood flow in human circulatory system 
in those vessels whose diameters are very small less 
than 500 pm, lpm=10'^m. 

The effects of peripheral layer, yield stress and 
other rheological Parameters on blood flow have been 
analysed. We have studied the problems on steady flow 
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of blood, effects of mild stenosis on blood flow through 
narrow vessels. Blood flow in narrrow tapered tube. 
Flow of Bingham fluid in capillaries. Pulsatile blood 
flow through narrow vessels. 

This thesis consist of seven chapters. Chapter I 
contains the general intrduction of the subject tviz. 
circulatory system, rheology of blood geometrical aspect 
of vessels etc.) and a brief account of the relevant 
literature. 

In chapter II, rheological behaviour of steady blood 
flow in narrow vessels have been studied. Two layers 
flow model, the central core layer and peripheral layer, 
satisfing herschel Bukley constitutive equation and 
casson equation have been considered for the analysis. 
Velocity profile and flow rate have been discussed with 
respect to the peripheral layer thickness, yield stress 
and parameter n. 

From the graph we find that as increases ( p 
remaining) the same) the flux rate decreases rapidly 
till Cp = .6 and till it has fallen to about 5 percent 
of Qq and then it rises again. Small changes in can 
make significant changes in Q. 

In Chapter III, a study of casson fluid model of 
blood flowing through a narrow tapered tube has been 
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done. Two layers flow model, the central core layer 
and peripheral layer satisfying casson constitutive 
equation have been studied for tapered narrow tube. 
Velocity profite and flow rate have been discussed with 
the help of tables and graphs. 


The variations of pressure gradient and shear 
stress at the wall have been computed for different 
tapered angles anddifferent axial distances, it is observed 
that the pressure gradient increases as the tapered 
angle and axial distanee increases. The variationof shear 
stressat the wall with flow rate for different tapered 
angles have been studied. We finch that increases 
with the increase in suspension concentration! and 
tapered angle. 


In chapter IV effect of mild stenosis on blood 
flow throw narrow vessels has been studied. In this 
chaper a non Newtonian Heschel - Bulkle^- fluid model 
has been' studieci in the presence of mild stenosts. 
Analytical expressions for volume flow rate, apparent 
fluidity and wall shear stress at the maxumufn height 
of stenosis have been obtained. ■ ' 


rate 


In this chapter we 
(Q) with stenosis 


study variation of volume flow 


height 






and yield stress 


(p) for different values of n (=0.5, 1.0, 1.5) 


cr 

Graphs show that value of Q decreases as p 

and p increases for fixed value of n. From results it 
is observed that the flow rate is decreased from 
Newtonian case by 25% for n = .5, 13% for n=l; 7% 

g 

for n = 1.5 for fixed and p. For n= .10, the value 

''0 

of Q are decreased by 30% from those for no stenosis 
case for fixed n and p. • 

In chapter V, a study of Bingham fluid model of 
blood flow in capillaries has been performed. In this 
chapter a two layer fluid model of blood with no slip 
velocity at the wall both satisfying Bingham constitutive 
equation is considered. In this chapter velocity field, 
volume flow rate and apparent fluidity have been found. 
Flow parameters have been explained. Results show that 
apparent fluidity fa, increases with d (plasma layer 
thickness) and decreases with (yield stress). It is 
also ‘observed that velocity field decreases fastly with 
r in the peripheral layer whereas in the core region 
velocity decreases slowly. 

In chapter VI a theoretical study of pulsatile 
blood flow through narrow vessels is performed. In 
this chapter we have analysed the pulsatile blood flows 
throw straight and long narrow rigid circular tube 


whose radius is constant. It is assumed that flow of 
blood is axially symmetric and in the z derection. In 
this chapter, the characteristics discussed are the 
pressure, shear stress at the wall, flow rate and 
viscosity. It is seen that the effect of frequency 
parameter a on velocity, shear stress, plug radius and 
flow radius is very small in the range 0<a<1- 

It is also found that values of Rp, and Q (t) 
vary greatly with respect to time specially for large 
values of reduced amplitude A. 

In Chapter VII, a theoretical model of steady blood 
floM/ in narrow vessel has been performed. Inthis chapter 
it has been assumed that fluid in core region satisfy 
casson equation and in mari\ginal layer region satisfy 
Newtonian equation. Apparent viscosity of blood has 
been determined as a function of yield stress, vessel 
diameter and peripheral layer (PPL) thickness. 

Here we have taken the data of Charme and 
Kurland and shown the variation of Pgpp with constant 
in a table, assuming that the casson vicosity is related 
by general equation - 


Where H is function of hematicrit and K is an 
expirimental constant. We can obtain Pappoidifferent 
concertration of Hematocrit. 



Where H is function of hematocrit 
ex pi ri mental constant, We can obtain 
concentration of Hematocrit. 



and K is an 
Qtdiffcren+' 
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CHAPTER I 

General Introduction 


1.1 Introduction 

Biomechanics is the field which has been developed recently 
by cross fertilization of the mechanics and biological sciences, so 
that both are being profited and are utilized for the benefit of 
mankind. 

. Biomechanics is the study where in the principles of fluid 
mechanics are used to understand the problems in biology. 

In the present study we concentrate to the study of blood 
rheology and blood flow in human circulatory system. 

1.2 Navier - Stokes Equations for the Flow of 
Viscous Incompressible Fluid 

Let u (x, y, z, t), v (x, y, z, t), w (x, y, z, t) and p (x, y, z, 
t) denote respectively the three velocity components and pressure 
at the point (x, y, z) at time t Jn a fluid with constant density 
p and viscosity coefficient // . Then the equation of continuity, 
which expresses the fact that the amount of fluid entering a unit 
volume per unit time is the same as the amount of the fluid 
leaving it per unit time, is given by - 


( 1 ) 


\ 


du dv ^ dw _ Q 
dx dy dz 

The equation of motion are obtained from Newtons second 
law of motion which states that the product of mass and acceleration 
of any fluid is equal to the resultant of all the external body 
forces acting on the element and to the surface forces acting on 
the fluid volume due to the action of the remaining fluid on the 
same element. The equations of motion, known as Navier. stokes 
equations for the flow of a Newtonian viscous incompressible 
fluid are - 
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( 2 ) 

(3) 

(4) 


If the external body forces x, y, z form a conservative system, 
there exists a potential function q such that 

da 


da da 

x=— , y=- 


dx 

dx 

dp 

& 


dy 


z=- 


dz ' 


— (Q + p), y--^ 
dx ^ ^ dy 


-d 

dy 


(Q + p), 


d 


dz 


{a + p). 


So that P is effectively replaced by P + q 


1.3 Reynold Number of Flows - 

In equations (2) to (4) the terms on the left hand sides 

represent the inertial forces (mass x acceleration) while the three 

terms on the right hand side of each equation represent respectively 
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the body forces, pressure forces and viscous forces. If U is a 
typical velocity and L is a typical length the inertial forces are of 
the order fU^ I L and the viscous forces are of the order 
The radio of these forces is of the order 
fu^/L fu^L^ fUL 


Re 


/ilU / L- 
h 


LuU 


UL 

V 


Where v = ~ is called the kinematic viscocity of the fluid. 
Now the diamensions of p and fUL are given by- 


Stress 

Strainrate 


fUL = ML-^LT-^L 

Thus Re is a dimensionless number. It is called Reynold's 
number, after asborn Reynold who in 1890 showed that the fully 
developed Poiseuille flow in a circular tube changes from stream 
line or laminar flow to turbulent flow when this number, based 
on the diameter of the tube, exceeds a critical value of about 
2000 . 

In straight pipes, turbulence occurs when Re exceeds 2000; 
in curved pipes it may not occur even for Re = 6000. 

When Reynold number is small, viscous forces dominate over 
inertial forces. 

1.4 Basic Concepts about Blood, Cardiovascular 
system and Blood Flows - 


Force, per. unit, area 
Velocity j lenght 

MLJ-^ / l3 


T-' 

ML-'r' 


= ML-T 




(5) 

( 6 ) 


Constitution of Blood - 

Blood consists of a suspension of cells in an aqueous solution 
called plasma which is composed of about 90 percent water and 
7% protein. There are about 5x10^ cells in a millilitre (Icc) of 
healthy human blood, of which, about 95 percent are red cells or 
erythrocytes whose main function is to transport oxygen from 
the lungs to all the cells of the body and the removal of 
carbondioxide formed by metabolic processes in the body to the 
lungs. About 45 percent of the blood volume in an average man 
is occupied by red cells. This fraction is known as the hematocrit, 
Of the remaining, white cells or leucocytes constitute about one 
sixth or 1 percent of the total, and these play a role in the 
resistance of the body to infection; platelets form 5 percent of 
the total, and they perform a functiuon related to blood clotting. 



Two types of forces act on a fluid element. One of them is 
body force and the other is surface force. The body force is 
proportional to the mass of the body on which it acts while the 




surface force acts on the boundary of the body and so it is 
proportional to the surface area. 

Suppose F is a surface force acting on an elementary surface 
area ds at the point P of surface S. Let and ¥2 resolved 
parts of F in the directions of tangent and norm^b at P. The 
normal force per unit area is called normal stress and is also 
called pressure. The tange«tWforce per unit area is called shearing 
stress. Hence F^ is a kind of shearing stress and ¥2 is a normal 
stress. 

Viscosity - 

Viscosity is that property of real fluid as a result of which 
they offer some resistance to shearing, i.e. sliding moment of 
one particle past or near another particle. Viscocity is also known 
as internal friction of fluid. All known fluids have this property 
in varying degrees. Viscosity of Glycerine and oil is large in 
compansion to viscosity of water or gases. 

1.5 Shear Tensor and Strain Rate Tensor- 

The shearing stress between any thin sheets of a fluid is 
defined as 

Force 

X = 

Area 

U 

now 

Vo 

Where yg is the distance between two planes s.t. one is at 
rest while the other is moving with uniform velocity U parallel to 
itself. 
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The fluid lies in between these two plates. 


Hence - h' 


WhereAiis the constant of proportionality and is defined as 


viscocity. 


by e. 


Sometimes y" is denoted by also it is denoted by ^ or 


is called velocity gradient, 


In the relation when n is constant it is clear that 

shear stress varies with velocity gradient. 

U. - \J ( 


0 No slib Cowdibovi , ^ , s 

( Si:a.f>on.e-r^) 

1.6 Non Viscous Fluid - 

A fluid is said to be non-viscous fluid if it is incapable of 
exerting shearing, stress. The following have the same meaning. 
Perfect fluid, Ideal fluid, Inviscous fluid. Non viscous fluid. 
Frictionless Fluid. In this case t = 0, so that 








1.7 Viscous Fluid - 

A fluid is said to be viscous if it exerts shearing stress. 
Thus in this case 

T 0 SO that j.1 0 


1.8 Elastic Bodies - 

Elastic Body is one for which velocity gradient vanishes. In 
this case 
du 


dy 


so that 


1-1 = CO 


/( du 


1.9 Newtonian Fluid - 

A fluid is said to be Newtonian if its viscosity does not 
change with the rate of deformation. In this case the equation 

f— 1 

\dyj 


T = jJ.| 

du 


is similar to the equation y=mx, where y=x , m=!.i, 


x='^. Hence Newtonian fluid i-s represented by straight line. 


1.10 Non Newtonian Fluid 

A fluid is said to be Non-Newtonian if its viscosity varies 
with the rate of deformation, aU are variables in the equation 
'du' 


T = p 


dy 



IdjLoi. Fiu-icL 


C SViear si^oi'n 


Therefore Non-Newtonian fluid is represented by curve. The 
main classes of Non-Newtonian fluids are Bingham plastics, Pseudo 
Plastics and Dialants. 


Non Newtonian fluids are those where stress is Non linearly 
proportional to strain (velocity gradient) 


1.11 ViscocIty of Blood - 

Blood is neither homogeneous nor Newtonian. Plasma in 
isolation may be considered Newtonian with a viscosity of about 
1.2 times that of water. For whole blood, we can measure effective 
visocity (apparent viscosity), and this is found to depend on shear 
rate. The constitutive equations proposed for whole blood are as 
follows : 

i) T = j-te" (Power law equation) 

ii) T = j.ie (Newtonian equation) 

iii) x = pe'’ + Xo(x>T o) (hershel - Bulkley equation) 

iv) (Casson equation) 

v) x = xy-4-pe where p is newtonian viscosity (Binghan 
equation) 



1.12 Cardiovascular System - 

The cardiovascular system consists of the following : 

1) The heart (which acts as a pump, whose elastic, muscular 
walls contract periodically making possible the pulsatile 
flow of blood. 

2) The distributory system (comprising arteries and 
arterioles for sending blood to the various organs of 
the body) 

3) The diffusing system (made up of fine capillaries which 
are in contact with the cells of the body) 

4) The collecting system of veins (which collects blood 
deplected of oxygen and full of products of metabolic 
processes of the system). 

The organs which supplement the function of the cordivascular 
system are (i) the lungs which provide a region of interphase 
transfer of O 2 to the blood and removal of CO 2 from it, and (ii) 
The kidney, liver, and spleen, which help in maintaining the chemical 
quality of blood under normal conditions and under conditions 
of extreme stress. 



Deoxygenated blood enters the right atrium (R A) from where it 
goes to the right ventricle (RV) when the heart contracts, the 
tricusped valve between the RA and RV closes and blood is pushed 
out to the lungs through the pulmonary artery (PA) which branches 
to the right and left lungs where (O2 is removed and blood is 
oxigenated. The blood returns from the lungs through the 
pulmonary vein (PV) to the left atrium (LA) and then it goes to 
the left ventricle. (LV) and from there, due to contraction of the 
heart, it enters the arota from which it travels to other arteries 
and the rest of the vascular system. 

1.13 Fahraous Lindqvist Effect (F-L effect) - 

In arteries, blood flows in two layers, a plasma layer near 
the walls consisting of only the plasma and almost no cells and 
a core layer consisting of red cells in plasma. 

1.14 Blood os a Transport AAedlum for the Body - 

Blood, the red fluid of the blood vessels with which we are 
all familiar, is the transport medium of the body. It is the medium 
by which all living tissues are related to their external environment; 
i.e. from the outside world. Blood is pumped to and fro between 
tissues and elementary canal to take up oxygen and nourishment 
respectively, and it helps in returning the products of oxidation 
and metabolism in the tissues to the outside world by the lungs, 
the skin and kidney. It is also the medium by which growth and 
repair substances are transported to the tissues and by which 
the various controlling glands of the body can distribute their 



Blood is some what viscid fluid. In man and in all other 
vertebrateanimals with the exception of two (the amphioxus and 
leptocephalus), it is red in colour. It consists of a continuous 
yellowish aqueous phase, the plasma, in which formed elements 
are suspended. The formed elements consists of red blood cells 
(erythrocytes), white blood cells (leukocytes) and Platelets 
(thrombocytes). Plasma is made up of water (92%) and contains 
traces of inorganic and organic salts. The inorganic and small 
organic molecules contribute little to plasma viscosity which is 
primarily dependent on the plasma proteins. If blood is allowed 
to clot and the solid material is removed, the remaining fluid is 
called serum. This has essentially the same composition as plasma 
except that fibrinogen and some of the clotting factors have 
been removed. About 7% by weight, of plasma is protein mainly 
albutnin, globulin and fibrinogen; which have molecular weight 
ranging from 44,000 to 10,00000. About half of the protein 
mass is albumin, the significance of plasma protein is its multi 
natures mentioned below : 

1- Plasma Proteins are responsible for osmotic pressure, 
the level of which is important for regulating water 
exchange between blood and tissues. 

2- They posses buffer properties and maintain the acid- 
base equilibrium of blood. 

3 - They produce a definite viscosity of plasma which is 
important in maintaining blood pressure. 

4- They promote stabilization of the blood by providing 
conditions that prevent sedimentation of erythrocytes. 
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5- They play an important role in coagulation. 

6- Plasma proteins are important factors in immunity. 

If blood to which an anticoagulant has been added, is poured 

into a rest tube and centrifuged, the corpuscles tend to settle at 
the bottom and the blood is divided into two layers, namely, a 
red lower layer consisting of the formed elements and a transparent 
colour-less or slightly yellowish upper layer consiting of plasma. 
The Leukocytes form of this white film between the erythrocytes 
and the plasma since their specific gravity is less than that of the 
erythrocytes. 

The red cells surface carries a net negative charge but in 
stationary blood, cells interact with each other to form aggregates. 
The aggregates commonly consist of 6 to 10 cells stocked face to 
face and such a cluster of cells in called a rouleaux. Secondary 
aggregation of rouleaux also accurs building up a complex three 
dimennsional structure. When blood is sheared, these rouleaux 
break upar>clat ^uffLuciently high shear rate the cell exists as an 
individual. Blood flows in the form of different laminar containing 
different types of cells. Considering these factors, the Rheological 
properties of blood might be expected to be rather complex. 

In normal blood cellular components others than red cells 
compose only about 3% of the total cell volume and usually have 
a small effect on blood viscosity. White cells may significantly 
inscrease blood viscosity in deseases where there is large increase 
in Leukocytes concentration. 

The red blood cells are about 97% of the total cell volume 
in the blood. Consequently, the removal of white blood cells and 
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platelets does not measurably modify the experimentally determined 
flow properties of these suspensions. The concentration of most 
red blood cells in suspension is generally reported as the suspension 
volume fraction occupied by the red blood cells, called hematocrit. 

The hematocrit is normally about 42-45% of volume. The 
erythroytes consist of a thin flexible unstrechable membrane with 
an interior filled with saturated hemoglobin solution (Viscosity 
6.0 centipose) The membrane is highly deformable if the change 
in the surface area is small, but becomes very much stiffer if the 
deformation produces a large change in the area of the cell 
membrane. 

The red blood cells have the shape of biconcave disc which 
can deform, however into a bullet shaped entity during passage 
through small capillaries. The cell has a large surface area relative 
to its volume. 

Two properties of blood are essential for the preservation 
of life. The first is that blood remains fluid in the blood vessels 
throughout its life and second is that it rapidly becomes solid 
when shed. The maintenance of fluidity is necessary for the 
circulation of blood, while the solidification of the shed blood 
provides indispensable defence against excessive bleeding from 
wounds. The coagulation of the blood is due to formation of a 
jelly by the deposition of protein material called fibrin and it is 
the formation of this body that the fundamental change in blood 
clotting occurs. 

To prevent this from occuring, various anticoagulents are 
added to the blood when it is drawn from the animal. 

The effects of an anticoagulcnt on rheological properties of 




blood are difficult to assertain because of the difficulty in 
making reliable rheological measurements in the short time available 
between the blood drawing and the blood clotting. The effect of 
anti coagulents on blood viscosity appears to be small. 

1.15 Human Circulatory System - 

In the human body circulatory system of blood consists of 
four functions Respiratory, Nutritive exetretory, Protective and 
Regulatory. The heart that provides the energy for the circulation 
has four compartments (ie left and right atrium and left and 
right ventricles) interconected to each other by one way valves. 
Blood coming from body tissues enters the right atrium. Through 
the venul caval contraction of right atrium forces blood pass the 
tricuspid vdlves into the right ventricle. From this point there 
begin two subdivisions the pulmonary and the systematic circulation. 
As the term shows, the former services the lung and latter the 
various systems of the body. The blood vessel have internal diameter 
in the range of 2.5 cm in aorta to about 4 microns in capillary. 

Blood is pumped from the heart into the aorta from where 
it goes to the circulatory system consisting of about 40 large 
arteries, 1600 main artery branches, 1800 terminal branches, 
4,00;00,000 arterioles, 1,20,00,00,000 capillaries, 8,00,00,000 
valves, 1,800 terminal veins, 1600 main venous branches, 40 
large venous branches, 40 large veins and then returns to the 
heart through the venule caval. 

It is estimated that the heart muscles it self consume about 
18% of the energy required to sustain life. only some energy 
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gose in to the mechanical work of pumping blood. The heart 
beats about 70 times in one minute in an average person at rest. 
Blood flow in smaller blood vessels approaches steady flow condition 
otherwise pulsatile. The pressure in the aorta rises rapidly to its 
maximum (systolic) value of about 120 mm/Hg. The requirement 
of the circulation is the supply of oxygen required for metabolic 
process. In an average man at rest, the 02required is 200 ml/ 
minute for under physical stresses, the need may rise to 5 litre/nrun 
.Iy\ an average man it is necessary for the heart to circulate 5 to 
6 litres of blood per minute or about the 5 times the metabole 
requirement of oxygen. The metablic requirements are the energy 
requiimients for biological functions. 

In fluid mechan-iGS the circulatory system is classified in to 
two parts v.iz.. macrocirculation and microcirculation. 

1.16 AAacrocirculation - 

The flow of blood in vessels of diameters greater than 500 
ym (such flow occcures in aorta) is called macrocirculation. The 
flow is characterized by high Reynolds number, defined as the 
ratio of inertial to viscous forces. Turbulent flow can occur in 
blood vessels for the Reynolds number greater than 2300. The 
governing equations for analysing of such flow conditions should 
Include inertial effects, effect of curvature of blood vessels, pulsatile 
flow and the distensiblity of the vessel wall. Blood can be taken 
as homogeneous and continuous fluid when vessel diameter is 
large in comparison to the dimensions of red blood cells. Pulsatile 
flow effects and Pressure wave propagations have been studied in 
macrocirculation by a number of research workers, Skalek and 
Taylor. The flow disturbances at bifurcations, bends etc. and 
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their effects on Pathological states are studied by Patel 


1.17 Microcirculation - 

when the diameter of blood vessel is less than 500 p. m [i.e 
arterioles, capillaries, venules] the circulation is 
microcirculatory.This is responsible for 80% pressure drop in 
circulatory system. In capillary bed the transfer of nutrients to 
and'removal of wastes from the living cells of the body is a part 
of microcirculation. The flow is characterized by verylow Reynolds 

Rn-r u/P ran not neelect the size of the red blood 


called 


1.18 Boundary Layer Flow, The Cone Flow and Fully 
Developed Flow of Fluid - 

when a fluid enters a tube from a large reservoir where 
the velocity is uniform and parallelto the axis of the tube, the 
velocity profile is a flat surface at the entry. 









Immediately after entry, the velocity near the surface is 
affected by the friction of the surface, but the velocity profile 
near the axis still remain flat. As the fluid moves further in the 
tube, the flat portions decrt^es, and at the section corresponding 
to A, the paraboloidal velocity profile for the fully developed 
flow is reached. The flow in the region OA is called the entry 
region (or inlet) flow and the flow beyond A (in region III) is 
called the fully developed flow. The length OA is called the entry 
length. The flow, in the entry length portion it self consistsof two 
parts. The flow in region I near the surface is called the boundary 



layer flow, the flow in region II is called the core flow or the 


uee 


In plug region, e = y 


£r T < xy where xy is yield stress 
£r V = constant 

plug flow exists where ever the shear stress does not exceed 
the yield stress. 

1.19 Rheology of Blood - 

Blood behave like a time dependent non Newtonian fluid 
and the basic rheological property of blood is its viscosity, the 
viscosity of hole blood I's about 4.0 centipose and of plasma is 
about 1-2 centipose at 37°C. The specific gravity of whole blood 
is 1.05 to 1.06 and that of plasma is about 1,03. Thus the red 
blood cells tend to sediment slowly in plasma. The increased viscosity 
relative to water is produced mainly by the presence of plasma 
proteens. Both, the molecular shape and concentration of the 
protein are important. Fibrinogens which has an elongated molecule 




has a marked influence; although forming less than 5% of the 
total plasma proteins. ■ ' ' 

.It is responsible for about 20% of the plasma viscosity 
elevations. 

1.20The Viscosity of Blood - 

It is an important factor in determining the local pressure 
variation through the cardio-vascular system which in turn 
influences the localflow rates through each section of the vascular 
network. The clinical importance of blood viscosity as a parameter 
lies in its sensitivity to small variations in composition. One can 
aften diagnose pathological states by determining a change in 
blood viscosity. 

There are several rheological parameters (e.g plasma, blood 
cells, hematocrit etc.) which can effect the blood viscosity, the 
viscosity of plasma increases with its protein concentration. But 
some proteins have different influences on plasma viscosity 
depending on their shape and size. The influence of Fibrinogen 
on plasma viscosity can be seen in the difference between plasma 
and serum viscosities. Serum usually has a viscosity which is 20% 
less than that of plasma. Many relationships have been suggested 
to express blood viscosity as a function of cell concentration, 
plasma viscosity and shear rate. When temperature is increased, 
the viscosity of blood and plasma is fallen. Measurements should 
be made at constant tamperature (37°C). 

Platelets and white cells, in general have little influences on 
blood viscosity, because they are present in very minute quantity 
as compared to red cells. ■ > d ' ■ ' 
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' V The deformation of red blood cells alloow the 

blood to remain fluid up to hematocrit of 98%; rigid cell without 
deformation will cease to flow at a rigid cell without deformation 
will cease to flow at a cell concentration of about 60%. The 
hematoerit also influence the deformation of the erythrocytes i.e 
raising of the cell concentration produces an increase in cell 
deformation and therefore fall in viscosity of blood occurs. 

In stationary blood, rouleax is formed which intracts to 
produce larger aggregates. At low flow rates the presence of these 
red cell structures strongly influence the viscosity of the blood. 
The siza of the rouleaux and aggregates progressively decreases 
as the. shear rate increases. In normal blood the disaggregation 
is probably complete at shear rate 50 sec.'l approximately. The 
viscosity of blood increases when shear rate falls, but it is 
unceration as to what happens when the ^hear rate actully falls 
to zero. The build up of a 3-D of increasing aggregates suggests 
that blood may show a yield stress. 

1.21 Blood as a Non-Newtonian Fluid - 

Blood, from fluid mechanics point of view can be thought 

as a suspension of erythrocytes in a Newtonian l|cjuid called plasma. 

In small vessels the dimensions of red blood cells are not negligible 

as compared to the size of tube. Therefore blood can not be 

considered as a homogeneous fluid in microcirculatory systesm 

(diameter less than 500 m,) It is a specific property of blood 

that exhibits a yield stress. Hence if applied shear stress is below 

a critical value (x <ty ) the response will be elastic and on removal 
of stress the shape of blood film is unaltered, othesrwise (T>Ty) 


-- . I.'.™. _*L.i__: 




flow takes place and blood behaves as a non Newtonian fluid. 

It is found that suspended blood cells are responsible for 
the non Newtonian nature of blood rheology. At low shear rates 
the blood exhibits yield stress and behaves like casson fluid 
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■[ x 2 -ty 2 + j.t 2 g 2 ] constitutive equations suggested by Herchel and 
Bulkley (T = Ty-rLie) has also been used to describe the shear rate 
dependence of blood. 

Blood rheology can be helpful in diagonosis of some blood 
problems. 


1.22 Survey of literature - 

The German form of the word (rheologic ^ and the description 
of srhall viscometer as a microrheometer are found in literature 
since early days. But in its modern sense, the term rhchology was 
coined by prof. E.C. Bingilam and formally adopted and defined 
at the foundation meeting of the American society of rheology in 
December 1929 in Washington, as " The science of deformation 
and flow of matter" 

However, the subject of hydrodynamics andaerodynamics are 
not included in rheologhy. 

The first mathematical paper of blood flow was given by 
Leon.-{^ard Euler in 1862. He developed one diamensional equations 
of invisicid flow of an incompressible fluid in an elasitic tube. 
Blood flow in captUaries is of great interest to Phisiologists involved 
in microvascular research. The pressure drop relation in 
microcirculation was oobtained by Posiseuille in 1940. The most 
remarkable result on blood flow was of Fahraeus and 
Lindguist, showing that the apporent viscosity (effective viscosity) 








of blood decreases as tube diameter decreases from 500 t-im. 
This effect has been confirmed by other investigators (Coklet 
and Dintenfass)0. Taylor studied dispersion process in New-dronian 
fluid flows- through a circular tube. 

,1.23 Geometrical aspect of Vessels - 

A striking characteristic of the circulatory system is its 
geometrical comple.,X4.ty. Blood must flow through many bends, 
bifurcations, stenosis and tapering during its journey, Most of 
the progress in understanding the mechanics of the circulations 
is based on the investegations of the flow in straight uniform 
tubes . 

It is only within the last two decades that the effects of 
geometric transitions on blood flow have begun to be explained. 

1.24 Stenosis - 

Stennosis causes the narrowing of the blood vessels due to 
the development of abnormal tissues and gives way to serious 
circulatory disorderby reducing the blood supp'ly. Hemodynamic 
charaterestics may be changed due to this undestTable growth 
which could be injurious to normal health. In recent years many 
workers have investigated the flow charcacteristies of blood through 
artery in presence of mild stenosis. 

Young gave a theoritical analysis of the effect of time 
dependent stenosis on flow characteristises of blood. Recently 
Sinha and Singh studied the effect of stenosis on blood flow 
through the couple stress.. 


1.25 Tapering - 

The mammalian arterial consists of the main aortic tube 
which continues into the artelaries and a number of side branches. 
Transformations and propagations of pressure and flow waves 
depend on the destribution of the. characteristic impedence. The 
two contributing factors are the decrease of, the cross sectional 
area and she peripheral increase of the valve stiffness. Many 
blood vessels have the characteristic of taper cylinders rather 
than straight cylinders. Block gave the idea that vessels which 
carry blood towards the tissues should be considered as long, 
slowly tapering cones rather than cylinders. The relations for 
shear stress and their variations with suspension concentration 
and tapered angle, have been discussed in this thesis. 









Table 2 

Mean Pressure (mm Hg)in Blood Vessels 






Protein 

Table 3 

Concentration in human 

Plasma 


Cone. 

Molecular 

Mol. Sizes 


(g/UKe^ 

Weight (ml.) 

■ (MJm) 

Protein 

75 



Albumin 

45 

^^,000 

15 X 4 

Globubin 

27 

1,50,000 

24 X 4.5 

Fibrinogen 

3 

3,40,000 

70 X 4 
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Chapter H 


Rheological Behaviour of Steady Blood Flow 

in Narrow Vessels 


2.1 Introduction 


From phisiological aspect blood flow throw narrow circular 
tubes is of great interest in microcirculatory system. From 
mathematical point of view blood is considered as a suspension 
of erythrocytes (red blood cells) in a Newtonian fluid called plasma 
when the diameter of red blood cells are no longer negligible as 
compared to the tube size. The apparent (effective) viscocity of 
blood depends on many factors such as, plasma viscocity, hematocrit. 
Size of vessel, shear rate, velocity of flow, rigidity etc. Further, 
the apparent bloo4viscocity in capillaries is much lower than that 
in large vessels (Dintenfass-) It is a well Known fact that apparent 
viscocity of blood decreases as the tube radius decreases . This 






is known as the Fahraeus Lindguist effect. Many workers in the 
field have shown that the fluid dynamic properties of blood flow 
like velocity profile, velocity gradient, shear stress at the wall, 
apparent viscocity of blood, deformability of red cells could play 
a vital role in the. fundamental understanding, diagnosis a.nd 
treatment of many deseases (cardio- vascular, diabetis etc). Many 
mathematical models have been provided for blood by a number 
of research workers, e.g. P.Chaturani, Nand Lai Singh, Charm et, 
al. 

Bugliarello and Sevilla have proposed two fluid models, in 
which either of both layers (Peripheral plasma layer and core 
hematocrit layer)l are of Newtonian fluid with different stress 
and vicosities. 

It has been found that for blood flowing in narrow tubes( 
500 pm- 20 pm), a cell free plasma layer occurs near the wall, 
whose diamension effects the blood behaviour to a large extent. 

Scott Blaier and Spanner reported that in normal range of 
shear rates (except at very high or a very low shear rate). There 
is no difference between Casson plots and Herschel Bulkley plots 
of experimental data and in some cases latter model provides 
better results. 

Many authors have the view that the slip at the solid surface 
has no evidence. Therefore we have neglected the slip velocity at 
the wall in this chapter. 

In our analysis we have considered the two layer flow model, 
constitutive equations of motion and the expression for velocity 
fields, flow rate and variation of flux with ~Ciobtained. 
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In the present analysis we have considered as incompressible 
viscou-sfluid model of steady blood flow throw narrow vessels 
under a uniform pressure gradient. It is assumed that the radius 
of vessel is constant. - 

Governing equations in cylindrical polar cordinate system 
(r^ where the axis of symmetry is taken as the axis of Z. 

There are, in general, three components of velocity, namely, Vr 
along the radius vector perpendicular to the axis, Vo perpendicular 
to the axis and the radius vertor, and Vz parrallel to the axis of 
Z, for the axe-symmetric case, we take Vo = 0 and we also take Vr, 
Vz and p to be independent of (gr; In this case the equation of 
continuity and the equations of motion are given by 


Hagen Poiseuille Flow 


Velocity Profile for Poisevllle Flow 






We consider steady flow when there is only one velocity 
component parallel to the axis so that Vr=o, V'gr^O and Vz=v. 
Then the equation of continuity gives 



The equations of motion (2) and (3) give 


must be a constant. Let us denote 


From equation (5) 


this constant pressure gradient by P. then 


Integrating this equation 


But velocity on the axis (i.e. at r=0) must be finite 
A=0, Then at r = a, V = o 


The total flow across any section i.e. total volume of fluid 
crossing any section per unit time is given by 


v2Kr dr 




Tta^P 
. 8|-i 


(on integration) 


.-.Q^ 


Tta^P 

8j.i 


( 9 ) 


2.3 Steady Non-Newtonian Fluid Flows in Circular 
Tubes 

Basic Equations for Fluid Flow 

' we consider the laminar flow of a non-newtonian fluid in a 
circular tube under a constant pressure gradient. Let the control 
volume be bounded by two coaxial cylinders of radii r and redr 
and let it be of unit length. Due to pressure gradient, there is a 
forward force px2 Hr dr on it. Let the stress be < 0 ) at a distance 
r from the axis. Then the force on the inner cylindrical surface 
is 2nn, and the force on the outer cylindrical surface is 

2n(rc) + ^Hd)r 

. Balancing the forces in the axial direction on the control 

volume, we get 


2Ti~{rx) = 2T:rP 

djr 




or = 


Integrating (1), we obtain 


rx = ^P^ + A or ^ 


1 


D 


( 1 ) 


( 2 ) 


Since the stress is finite on the axis (i.e. at r-O) We have 


A=o, D = 0, X= 2 *^*^ 


C3) 


The velocity v is parallel to the axis. It is a function of r 



only ard is expected to decrease from a maximum on the axis to 
zero on the surface so that the only non zero compohent of 

strain rate is 


e = 


fcJv 
'I dr 


For a non-Newtonian fluid 


xc=f(e) 

From (3) “ (5) 


1 / dv 

Ipr f 

2 V dr 


(4) 


(5) 


( 6 ) 


Integrating (6) subject to the condition that v=o when r-R, 
we get V as a function of r, Then we can obtain the flux Q by 


using 


R Cj') 

Q = ^2%rvdv ^ 

Integrating right-hand side of ( 7 ) by parts we get 

Q = 


Vvl ^ 

' J 2 dr 


' o J 2 dr 
0 0 

Since v = o at r=R we get 

Q = 


( 8 ) 


(9) 


;.4 Flow of Power-low Fluid in Circular Tube 


Here x = j-te 


dr 


"ml" 

r 2 i-i 


Integrating (10) we obtain 


( 10 ) 









In this case, we have e-0 when t<t„ and there is a cot 
negion which flows as a piug. Let the radius of this plug region 
be rp. At the surface of the plug, the stress isTo so that, considering 

the forces on the plug, we get . 

PxTcrp^ = 

2to (1) 

or — 

In the non-core region, x>Xo and 


X = i-ie" 4- Xc 


Integrating (4-) we obtain 


^ Jl^i 1 " {R - rp)r' - (r - rp) 
n + 1 v2iJ L 


If r=rp then v=vp (the velocity of the plug flow) so that 


” P \ (R-mV" 

Egu (1) determines the radius of the plug and then using 
this value of the plug. (6) determines the velocity of the plug and 
(5) determines the velocity in the non-core region. 
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Results & Discussions - 

This gives the relative change in Q with To, Figure. 2.4 
illustrates the vaiations of f(Cp) with Cp for various values of n. 
The figures shows tho± 


Tirp^vp + jZntvdr 


" ( Ji(. ^ -p) - ^ 

(n + 1 2^ ^ ' ■ 1 + 3 -Us 

^ ^ n n J 


n (PypT^ 

(nul)l2).Lj 


C/(1-C,)”*’+(1+C,){1-C,)=* 

2_ + 3 - + 2 . 

n n 


n r PV oU 


(3n +1) v2p 


R" Hj(Cp) (say) 


Where + = -^ = |i,f(0) = 1 


IfQodenote the flux when there is no plug flow (j£e where t -0 
Cp=0) we get 









2.6 Flow of Casson Fluid in Circular Tube 


Here 


1 1 
^2 


+ To2(t > Tq) 

So that for the noncore region 


( 1 ) 


-dv 

dr 


e = 


1 Y2 

2^ 


1 \ 2 
12^ 


1 / 

JJ-,; 


or 

dv 


P 

M 


/ 1 

~rp‘ 




YP 

2 (.1 


(2a/w- 


r-AP 


Integrating (2) we obtain 


2 p 


3 


+ rpP 


( 2 ) 


(3) 


so that the plug velocity is given by 


=I£[1r 2 - p_i;^p2 -Irp 

2pU ^ S'/p 6 


1 PR^ ^ 8^1-^ 

I p 3 p e p ; 


4 Ll 
PR" 


4p 


g(cp)(say) 


v,. 


Then, 




P(cp) 


(P) 


Results a?id Discussions 


Fig.2.>a- shows the variation of g (cp) with cp. This shows 
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that as Tb increases CM- remaining the same). The plug velocity 
or the maximum velocity of flow decreases rapidly till cp reaches 
.6 when the velocity is reduced to about 6 percent of the value 
and then it rises slightly. For blood, samll changes in Tt* lead to 
significant changes in maximum velocity. 

The flux Q is given by 


ToThat 


Results and Discussions 


Figure. . the graph of h (Cp) against Cp. 

It shows that, as To increases (/4. remaining the same), the 
flux"deGreases rapidly till Cn= .6 and till it has fallen to about 5 


‘ V'; 


J '' J' ' 




percent of Qo and then it rises again For blood, small changes in 
Tb can make significant changes in Q 

The Casson fluid flows in the tube only if iKp<Ft ii-e if 
I ■ 2”^o P R 



Vari’a-HeYi of 

a., a 
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Chapter Mi 


A Study of Cason Fluid AAodel of Blood 
flowing through a Narrow Tapered Tube 


3.1 Introduction - 

The behaviour of blood flow is influenced by the suspension 
of cells in an aguas solution of organic and inorganic substances 
called plasma. There have been many attempts to explain the 
anamalous behaviour of blood by proposing defferent theoretical 
models, e.g. Newtonian and non Newtonian ( microplar casson, 
Herschel- Bulkly, Power lawer etc.) 

lida and Murata (64) have studied pulsatile blood flow through 
small vessels by assuming Hershel Bulkly fluid model of blood. 

Ariman et al. (3), Bugliarello and sevillaj Chaturain and 
Llpadhya and many other workers have investigated the flow of 
blof;d through straight circular tubes. Since the vessels bifurcate 


36 


to distance. Due to the reduction in the vessels diameter at each 
bifurcation and angle of taper (which is small) thier influence 
can not be neglected 

Block (8) has suggested that all the vessels which carry 
blood towards the. tissues should be considered as long, slowly 
tapering cones rather than cylinders, Chaturani and prahlad have 
investigated a steady laminar flow of blood in a uniform tapered 
tube by considering blood as polar fluid. 

Many research workers like Dwivedi et al have studied steady 
blood flow in tapered tube by assuming different models of blood. 

Some research workers analysed both theoretically and 
experimentally the aspects of blood rheology in tapered' tubes 
and obtained pressureflow relation 

Block and oka studied pressure loss of blood flow through 

tapered tubes. 

Hence in order to have a better understanding of the complex 
nature of blood flow, one would not only take in to account the 
^properties of blood but also the geometrixal shape of its vessels. 
In present chapter, we have studied the flow of blood through 
uniform tapered tube assuming that it obeys the casson fluid 
.Vnodel. The expressions for pressure gradient and wall shear 
stress have been obtained. Variations of pressure gradient and 
shear stress at the wall have been studied for constant suspension 
concentration with respect to tapered angle over the flow rate 
range 0.01 to 0.1 cc{ sec Numerical values of pressure gradient 
for Bingham, Power law and Newtonian fluids have been obtained 
for fixed value of suspension concentration. 



3.2 Mathematical Model 


&e.ohief 


Let US consider a steady laminar flow of an incompressible 
Non-Newtonian, casson fluid in a uniformly tapered vessel is shown 
in above fig 3-‘l-.The geometry is most conviniently represented 
by the-cylindrical Coordinates: The radius R (z) of the tapered 
tube is given by 

R .(z)= Rq ■ ^ 6- (1) 

where Rq is the radius at Z=0, G- is the tapered angle and 
Z is along the ay:\s of the tapered tube.The problem has been 

studied under the following assumptions. 

The motion has an axial symmety, 

^ . No body forces act in the fluid. 

(3)" the motion is so slow that the inertia term can be 
neglected. 





where P, and V are respectively pressure , shear stress 
normal to r in z direction and the axial velocity. 

The casson constitutive equation is given as 


where )/- — > 'Vy and i-i represent shear strain, yield stress 
dr ■ 

coefficient^viscocity rep actively. 

The oppropriate boundary conditionsare given by 

V =0 at T" = R (z) 


and is finite are r =o ' 

Where Ry is the plug radius and Vy is the plug velocity 


3.4 The AAathematical Analysis 

In the analysis to follow, it is assumed that the pressure 

gradient is a function of axial coordinate only. Then with the 
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help of integration of equ. (2) with boudary conditions, yields 
further, 

r dp 


z dz 


From equ. (5) Ex (10) we have 


( 10 ) 


1 - - 

Tq 2 + 


2L 


Or y 
hence 

dv _ ^ 
dr ri 


VilJ l2L ) 


xw 


1 

^2 




-Xo2 


( 11 ) 


Here L is the length of the tapered tube. Integration of 
(II) with the help of boundary conditions, we get 


t (ii ™ , , 

v = ^R{z) 

2ii 


Where P 


R\^) ^ 


^ X ^ 


i ^ 

i_JfL 
I R\2). 


+ 2p 


^ 7“^ 

V RqJ 


( 12 ) 


On substituting r=Pf?(z) . the plug velocity is given by 

8 J 


Vy=^R{Z) 


1-;p"^+2|3-l|3^ 


3" ‘ 3 

VoLulretnc Flow Rate 

The volume flow rate Q is given by 


(13) 


Q = Qy + Q 


(14) 



where Q- 


nV„R. 


Now substituting for V and Vy from gqu-(12) £r (13) into 
equations (15) £r (16) respectively, we get on neglecting the 
higher powers 


Hence 


The pressure gradient is obtained from equ. (10) £r (19) 
with the help of boundary condition (6) as 


The shear stress at the wall is given by equ. (19) as 

2riQ -j , ® «2 R ^ /'> 1 

Finally equ's (20) Sr (21) yield 
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Results and discussions - 

It may be noted from equations (20) and (21) that the 
pressure gradient and the shear stress at the wall are not constant 
along the tube axis. 

Infact, they increase with the decrease in the radius of the 
tapered tube. From equations (22) it is observed that the shear 
stress at the wall depends on the pressure gradient. Dwivedi et 
al have proposed the micropolar fluid madel for blood flow through , 
a small taperd tube, where they have assumed the existence of a 
constant pressure gradient throughout the analysis. It appears 
to be incorrect for the reasons prescribed above. 

The variations of pressure gradient and shear stress at the 
wall are computed from egus (19) to (21) with flow rate over the 
range .01 to .1 cc/sec. for different taperd angle amd defferent 
axial destance. It is observed from fig3'‘145'2that the pressure 
gradient increases as the tapered angle and axial distance increases, 
flow rate being taken over the range 0.01 to 0.1 cc/sec Figure 
.....SiS'. show^the variation of shear stress at the wall with flow 
rate for defferent tapered angles. We find that ijvxreases 

with the increase in Suspension concentration and tapered angle. 

■ tWj'^also an increasing function of axial distance.. Therefori* 
for known flow rate, the shear stress can be computed at any 
point of the tapered tube. 

■ , , dp . 

Tables. .3.f... ( a, b)show that for Bingham model — increases 
with the increase in tapered angle and axial distance for a fixed 
suspension concentration 40% however, these values are 
comparatively less than those for casson model. 


Few Newtonian fluid, through ^ has the same trend as for 

uZ 

casson. its values smaller than those of casson model. It is hoped 
that the present model will provide a better understanding of 
vascular fluid mechanics. 


Table 3.1 

Bingham Fluid (n = l) 
dp 

(a) Variation of — with Q for defferent A 

dz 

H = 40% = .01 cm Z = .10 cm 



dp .. 
& 

nes/cm^) 

X 10 ^ 



Q(cc/sec) 

0=1.0° 

1.2° 

1.4° 

1.6° 

1.8° 

0.02 

0.537 

0.633 

0.755 

0.901 

1.097 

0.04 

1.096 

1.295 

1.540 

1.840 

2.283 

0.06 

1.634 

1.930 

2.295 

2.793 

3.33 6 

0.08 

2.193 

2.590 

3.080 

3.682 

4.477 

,0.10 

2.736 

3.233 

3.844 

4.595 

5.588 


References ^ 

1] lida and Murata, 1980: Biorheology 17, P.377 

2] Arimanetal, 1974: J Appl. Mesch. 41 PI. 

3] Block, E.H, 1962: Amer .J. Anat. 110, P.125 

4] Dwivedi et al, 1982: Ind. J. Technol. 20, P. 295 

5] Oka, s, 1964 : J phys. Soc. Japan 19, P 1481 


44 


dp 

(b) Variation of — with Z for different Q 


H= 40%, Ro = O.OlCm 0= 1.4 



dp 

dz 

(dynes/cm^) X 10 ^ 



Q 

Z = 0 

0.05 

0.10 

0.15 

0.20 

.02 

0.243 

0.414 

0.755 

1.571 

3.559 

.04 

0.497 

0.845 

1.541 

3.085 

7.262 

,06 

0.741 

1.260 

2.297 

4.998 

10.792 

.08 

0.994 

1.691 

3.083 

6.171 

14.524 

.10 

1.241 

2.110 

3.848 

7.701 

18.127 


Table 3.2 
Power Law Fluid 

(a) Variation of — with R for different 0° 


H=40%, Rq = O.Olcm Z=.01cm n= .75 


*' 

t .«» 

Q 

0=1.0° 

1.2° 

1.4° 

1.6° 

1.8° 

0. 

02 

2.062 

2.377 

2.726 

3.216 

3.775 

0. 

04 

3.107 

5.893 

6.760 

7.968 

9.348 

0. 

06 

6.899 

7.961 

9.132 

10.764 

12.128 

0. 

08 

8.556 

9.874 

11.326 

13.350 

15.662^ 

0. 

10 

10.124 

11.684 

13.401 

15.797 

18.532 

(b) 

dp . 

Variation of — with Z 
dz 

for different Q 



Q 

z=o 

0.05 

0.10 

0.15 

0.20 

0 

.02 

1.104 

1,706 

2.740 

4.429 

97.736 

0 

.04 

3.736 

4.218 

6.794 

10.944 

24.168 

0 

.06 

3.696 

5.698 

9.178 

14.784 

32.648 

0 

.08 

4.584 

7.068 

11.383 

19.336 

40.49 2 

0 

.10 

5.424 

8.362 

13.469 

21.696 

47.912 


Table 3.3 
Newtonian Fluid 


(a) 


dp 


Variation of — with Q for different 0 
8z 


H = 40%, 


Ro=.01 cm 


Z=.10cm 


Q 

0=1° 


1.2° 

1.4° 

1.6° 

1.8° 

0.02 

0.418 


0.496 

0.589 

0.703 

0.853 

0.04 

0.857 


1.009 

1.203 

1.428 

1.744 

0.06 

1.269 


1.504 

1.793 

2.128 

2.594 

0.08 

1.703 


2.019 

2.407 

2.850 

3.488 

0.10 

» 

2.125 


2.520 

3.004 

3.564 

4.313 

(b) Variation of 

dp 

dz 

for different values 

of Q and 

Z 

Q 

Z=0 


0.05 

0.10 

0.15 

0.20 

0.02 

0.23 2 


0.312 

0.411 

0.512 

0.523 

0.04 

0.435 


0.488 

0.512 

0.612 

0.711 

0.06 

0.812 


0.912 

1.008 

1.111 

1.212 

0.08 

1.805 


1.015 

1.189 

1.212 

1.314 

0.10. 

2.004 


2.810 

2.981 

3.005 

3.502 
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Table 3.4 
Casson Fluid 

dp 

(a) Variation of — with Q for different 0° 



4-0%, 

Ro =0.01 

cm 

Z =1.4 


Q 

00=10 

1.2° 

1.4° 

1.6° 

1.8° 

0.02 

3.001 

3.278 

3.415 

3.912 

4.775 

0.04 

4.201 

6.022 

6.819 

7.112 

8.992 

’’ •■■’0.06 

- 7. lbs 

9.002 

11.178 

■• ■13.5 2 7 

15.819 

0.08 

9.552 

11.553 

13.112 

15.778 

17.111 

0.10 

10.812 

12.212 

14.013 

16.123 

18.412 

dp 

(b) Variation of — 

with Z for different Q 


P = 

0.01, 

Ro=-01> 

• 

0 =1.2°. 

^=0.1 

Q 

Z = o 

0.05 

0.10 

0.15 

0.20 

0.02 

0.280 

0.303 

0.338 

0.382 

0.44 

0.04 

0.560 

0.607 

0.676 

0.765 

0.88 

0.06 

0.840 

0.910 

1.015 

1.148 

1.3 2 


Table 3 


dp 

(c) Variation of — with x (z) 



Ro-.Ol, 

z=.o 




o 

rH 

11 

O 

O 

1 . 2 ° 

1 . 4 ° 

1 . 6 ° 

0.02 

4.848 

5.063 

5.263 

5.556 

0.04 

9.696 

10.127 

10.526 

11.111 

0.06 

14.545 

15.189 

15.789 

16.667 

0.08 

19.394 

20.253 

21.052 

22.22 

0.10 

24.292 

25.316 

26.316 

27.778 



Chapter IV 


Effect of Mild - Stenosis on Bio 
Through Narrow Vessels 



A stenosis causes the narrowing of the bloodvessels due to 
the developement of abnormal tissue and it gives way to serious 
circulatory disorder by reducing or occluding the blood supply. 

It is well known that abnormal growth in the lumen of the arterial 
wall developeat various locations of the cardio-vascular system 
under diseased condition. 

Arteriosclerosis or stenosis, as it is called, is one of the 
most widespread arterial diseases. The actual cause of stenosis is 
not well known but its effects on the cardiovascular system can 
be understood by studying the flow behaviour and flow 
characteristics in its vicinity. Several theoretical and mathematicillL 
attempts have been made to study the blood flow characterists 
due to the presence of a stenosis in the lumen of a blood vessel. 

Assuming blood'' to be a newtonian fluid many authors like 
chakravarty and chawdary [22], Doffin and Chagneaw [41] Liou 
Ital [79], Mac Donald [82], Mann etal. [83], Mehrotra andjayaramah 
[88], Misra and chakravaty [92], Morgan and young [93], pollard 
[112], young [176] and young etal. [178] studied . the blood 
behaviour in stenosised tubfi. - 

Mehrotra and Jayaraman [88] have stucLled steady laminar 
flow of blood in an elliptic stenotic tube by assuming blood as 
Newtonian fluid and analytical expressions for shear stress 
distrubution and impedance are obtained. Misra and Chakravarty 
[92] studied blood floo) through an arterial segment having a 
stenosis and the artery is modelled as an initially stressed 
erthotropic elastic tube. They have assumed that blood is a 




nev/'i^tenian fluid In small vessels and ^'at low shear rates the 
flow behaviour isNlot Newtonian, and it has been seen that the 
velocity profiles are not parabolic under these flow conditions. . 
Bugliarello et al (16), Randetal (117) and many others have 
consideraed the non NewtoniCOi behaviours of blood in their 
studies. Recently sirw^aand Singh [130], Srivastava [138] studied 
the effect of stenosis on blood flow, through the couple stress 
fluid model of blood. 

Bitoun and Bellet [7], Rabinovets et. al, Slazback etal. 

[134], and Tandon et al [150]. have studed pulsatile flow of blood 
through rigid and pulsating stendsjes. i. ' ' ■ ■ 

, • ' ■ ' ■, Irhe geometry, of constriction chosen was defined 

by a Gaussion normal destri.bution curve, young and Tsai [177] 
studied some of the flow characteristics in models of arterial 
stenosis under steady and pulsatile flow conditions. These 
experiments have yielded, in the case of steady flow, a description 
of the extent of separated flow regions, a measure of pressure 
losses of constriction. Several other workers studied the 
characteristics of blood flow .due to presence of stenosis in the 
lumen of the artery [ Abdellah and Hwang] [1], Chakravarty [20], 
Devanathan and Parvathamrna [37], Nakameera and Sawada [95], 
Pralhad and schultz [114] sud and sekhon [141] and.yamaguchi 
and hanai [147] by assuming different models of blood, Further, 
Shukla et al [127] studied the effect of Peripheral layer viscocity 
of blood on the flow characteristics with mild stenos^Shukla et 
al[128] investigated the effects of stenosis on peripheral resistance 
and wall shear stress in an artery by considering the blood as 
power law and casson's fluid.. But it has been, reported by scdtt 
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- Blair and Spanner [148] that the hershel - bulkey fluid model 
of blood has an edye over other fluid models. Therefore studying 
the effect of an axially symmetric mild - stenosis on the flow of 
blood in narrew vessels, we, assume that blood obeys hershel - 
Bulkley constitutive equation of motion. The volume flow rate, 
apparent fluidity of blood a.nd shearing stress at the wall in the 
stenosized tube are discussed with respect to different field 
parameters. 


4.2 The Mathematical Models - 


In the present chapter we analyse a Herschel- 
Bulkey fluid model of blood flowing through narrow 
vessels with mild stenosis, cylindrical coordinates system 
( r, z, 0 ) has been considered with z - axis coinciding 
with the central line of the vessel. The geometry of 
the stenosis as in fig..ff.-.l.... is described as 


R. 


2R. 


1 + cos^[z-d-^ 


d < z < 4 + d 


( 1 ) 


= 1, otherwise 


We 


shall assume further that 




(a) (b) ^-0(1). (O Kf r^'' C 2 ) 

0 L.\mJ V y 

where Rq is the radius of unobstructed tube, R 
the radius of obstructed' tube, Lq the length of the 
stenosis, d the location of stenosis, the maximum 
height of stenostic growth and Re is the Reynolds 
number of the flow. 


4.3 The Governing Equatlon- 

The equation of motion and continuity for steady 
incompressible laminar flow under no body forces are 
given by 


dZ rdr 




Where P is the pressure. 


the shear stress 


normal to r in 2 direction, v the axis velocity. 

The constitutive equation in one-diamensional form 
of Hershel-Bulkey fluid is^^^^^^ 






f i, ' -ivi ^ ■’ if v? <;4 s:! 'ii' ' "T'm 






T = }.iy'’ +t:o;x>T(. 


y = 0 :t < Tq 


( 6 ) 


Where x represents shear stress, Xq the yield 

' dv' 


stress, pthe coefficient of viscocity and Y \dry 
strain rate. 


the 


4.4 The Boundary Conditions - 

The appropriate boundary condition are given by 


v=0, at r=R 


X is finite at r = 0 


(7) 


( 8 ) 


Where .x^is the shear stress at, the stenotft • wall 






•t - ■- 




4.5 The AAdthematicol Analysis - 

Integrating equation (3) we get 


dp r ^ 
rx= — — + c 
^ dz 2 


(9) 


Where C is constant of integration 

Using the boundry condition (8) in equ. (9) we 


1 


get - 


( 10 ) 


AP 

" 2L ^ 

Where aP 's the pressure drop over a lenght L 
of the tube. 

From equation (6) with boundary condition (7) 
we get - 




n f 


(n + 1) V p 


n+1 


n+1 




( 11 ) 


p = ^ 

where _ 

Plug flow exist whenever the shear stress does 
not exceed yield stress and the plug velocity Vy can 
be obtained by putting r=PR in equ (11) 


n ■ X. 


1 


n+1 


R(1- P) " 


( 12 ) 


(/7 + l)l. P 

From equ. (11) we can ^obtain velocity equation 
for Bingham and power law fluids by putting r=l and 
p = 0,' respectively. 

4.6 The Volume Flow Rote - 

The volume flow rate Q is given by 
•Q - Qy + Q’ ■ . , 

Where 



where 


(3/7 + 1) V f-i 


At the maximum /height of stenoses the volume 
flow rate is given by 


Qy \27iV/dr = nV^R^ 

0 

R 

Q> _ ^2nvrdr 


( 14 ) 


(15) 


Using equation (11) £r (12) in (14 Sc (15) 


n+1 


Ro^P^(1-pF 


(16) 




n + 1 V i-i 7 


Rn 




2il±l 

n 


2n^ 

3n+lf d\ 

(2/7 + 1)(3p + 1) 

^ UJ 


(17) 


Hence for equ. (12), (16) &c (17) we have 


nit . "o 3 1 

Q- 


3/7+1 

1 rx rp 


/?(2;7 + 1)^JVR0; (18) 

The non-diamensional flow rate Q is given by 


. 3p + 1 

1- p 

n(2/7 + 1) 
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Q = 


3n + 1 


n 


(2n + l)jV R, 


1 - 


3a 


( 20 ) 


^0^ 


VJ 

higher orders of 

4.7 Apparent Fluidity and Shear Stress 
at the Wall - 

From equ. (18) the apparent floidity <^a is given 


by 


= 


(3n + l)Q 


mzx^nR^ 


or 


^a = 




1 


'"rT 


( 21 ) 


higher powers of p are neglected. At the meximum 
height of stenosis i.e. at z = d + ^ we have 


(j)a = 


3n + 1 
2n + 1 


1-3n 


R. 


( 22 ) 


■o2 


From equation (18). the shear stress x, at the 
wall is given by ' 


(3n + l)Qp'’. 


rmRn 


^ 3n + 1(5 
1 + :— iP 


2n + 1 2 \Ro2 


^ r'^ 


-Zn 


(23) 


Now the wall shear stress at the maximum 


59 


height of stenosis i.e. at z = d + ^ is obtained, as 


(3n + l)QiJ,'’ 


rniRr, 


1 + 


3n + 1 
2n + 1 


P 


1 + 3n 


a 


V 


R, 


oJ 


(24) 


In non- dimensional form, we have 


Xw = 1 + 


3n + 1 
2n + 1 


r 


\ 


1 + 3n- 

V 


(25) 


Tw 


X = 


(3n-Hl)Q|i” 


miRn 


, x' being the Wall shear stress for 


-no stenosis case. 
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Results and Discussion - 

In the present chapter a non - newtonian Hersche! 
- Bulkley fluid model has been studied in the presence 
of mild stenosis Analytical expressions for volume flow 
rate/ apparent fluidity and wall shear stress at the 
maximum height of stenosis have been obtained Further, 
it is observed that in the earliar models including the 
two fluid models of Shukla et al [ 1 2^ ] the blood 

has earliar been assumed to be a Newtonian fluid or 
has been modelled as double layered Newtonian flu-id 
with djLfferent viscosities in the core and the peripheral 
region But erperimental results ( Cokelet [3+] and 
Goldsmith and skalak [ 53 ] clearly indicate that blood 
exhibits non - Newtonian effects and also indicate (for 
flow through small diameter tubes) the existence of a 
peripheral layer of Newtonian fluid and core of non - 
Newtonian fluid (RBC suspension) In our analysis it 
is assumed that blood is non Newtonian and flow is 
single fltiid model. The results are discussed with 
stenosis height^ yield stress of blood and non - 
Newtonian parmeter n ■ 

Figures and depict the variation 

of volume flow rate (Q) with stenosis height 
and yield stress (p) for different values of n (=0.5, 
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1.0, 1.5). The plot for n = 1.0 is for Bingham fluid. 

From figures it is clear that Q decreases as andB 

increases for fixed value of n. The values of Q are in 
increasing order with parameter n. From results it is 
observed that the flow rate is decreased from Newtonian 
case by 25%5for n =.5, 13%, for n = l; 7% for n = 1.5 

for fixed p and p. For p =0.10, the values of Q 

are decreased by 30% from those for no stenose case 
for fixed n and p. 

...OLinii- _ show^the variation 

of apparent fluidity (j)a vwith different 

value of PKfh(=. 5, 1.0 1.5). From graphs we observe that 
apparent fluidity decreases as stenoses height and yield 
stress increases for fixed n. The value decreases 
uniformly with parameter n. From numerical resultes, 

CJ ■■ ■■ 

we abtain that for = .10, the values of are 

decreased by 15% for n = 0.5, 30% for n = 1.0; 45% for 
n = 1.5 from those for no stenosis case. Further we 
observe that (j)a is decreased by 13% for n = 1.0, [5=0.10 
as compared to Newtonian, value. On comparing the 
results with Newtonian we find that at a fixed stenosis 
height, (t)a is greater than that for Newtonian value 
for n<l and smaller for n>l. 
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% , 


TaUes 


and .} depict the 


variation of wall shear stress with and p for 

different values ofP4’^ (=0.5, 1.0, 1.5). From figures, 

a 

it is seen that xw increases with respect to ^ and p 

for fixed n. The valus of are in increasing order 

with n. The results show that the value of t,„ is 

w 

a 

increased by 17% (approx) for n = 0.5, '^=0.10; 30% 

a <y 

for n = 1.0, ^=0.10; 45% for n = 1.5, =0.10 for no 

- stenosis case Also we obtain that xw increased by 
13% from Newtonian case for n = 1.0, p=0.10. 

Srivastava [ 13€ ] reported that for couple stress, 

fluid wall shear stress increased by 60% from no - 

stenosis case it is noted that the values of wall shear 

stress are greater than Newtonian value for n>l and 
smaller for n<l. Our results of xw 

conformity with the results of Srivastava [136 ] 

for couple stress fluid model. 

From equ 2a. and Zk we have 

calculated the numerical values of (t)a and Xw power 
law fluid (p=0) (Table ..1L. Be .2....). From these tables, 

. . . ■ CJ', ^ 

we abserve that decreases with increase in and 

g ' 

n, and xw increases with respect to ^ and n. 
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TA8LE4-3 

VARIATION OF WITH s-/R„ FOR DIFFERENT (3 


4a 


/ * 

0 

il 

CD- 

•15 

-20 

-25 

•05 

0-257 

0-177 

0.10 2 

0.0/5 

•io 

0-25 ^ 

0-1S3 

0.089 

0-01^ 

•15 

0-215 

0.148 

0.0?>1 

0- 013 

•2o 

0-19 4 

0-154 

0.073 

0- 012 
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T-38LE^-+ 

UfTH «-/R„ for different /S 
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CHAPTER V 


A Study of Bingham Fluid Model of 
Blood Flow in Capillaries 


Blood flow in capillaries of internal deameter equal 
to that of red cell is quite different from that in the 
large vessels. Since the size of RBC is not negligible 
as compared' to vessels diameter, it. is necessary to 
consider the flow of blood in capillaries as two phase 
non homogeneous flow. 

The apparent viscocity of blood depends on several 
factors such as plasma viscosity, hematocrit, size of 
the vessel, shear rate, rate of flow, regidity and 
deformability etc. Also the apparent viscosity in 





capillaries is much lower than in large vessels 
(Dentifass). Also it is well known that apparent viscosity 
of blood decreases as the tube radius decreases. 
(Fahraeus lindguist effect). Gupta and Seshadri, Lida 
and Murata, Charm et al have considered the two fluid 
models, in which both layers (Peripheral Plasma layer 
and core hematocrit layer) are of Newtonian fluids 
with different viscosities and yield stresses. 


5.1 The Mathematical Analysis 


P.4yu|>AAJ>LC»-i. g^p 


Jxe.^X9\/y. 




V 


The flow -ot blood is axially symmetric and is in 
the Z direction. 

The equation of motion and continuity for steady, 
incompressible, laminar viscous flow under zero body 
forces are - 


- dp 1 d 
dz r dr 




Where t is shear stress normal to r in z direction, 
P is the pressure and V is velocity. 

Let a two ' fluid model of blood with a central 
core region of radius Rc and a plasma layer of thickness 
J ^ (Rp satisfying Bingham plastic constitutive 

equation. 

( 4 ) 

X = Xj, + fXg, X > ty 

y =0> 

where ^ is Newtonian viscosity, x is shear stress, 
Xy is yield stress,, y is shear strain. 

Boundary conditions are 
Jp=Tc, Vp*Vc at r=Rc 


















..’.A- 1 '*' 


Vp=0 at r=Rp 
y = 0, = Vc=Vy at r=Ry 


(5) 



■ 


■r is finite at r 


0 


Where Ry = Plug flow radius, Vy = Plug flow 
velocity, suffixes P and C denote the value of the 
corresponding parameter in peripheral layer and core 
layer, respectively. 


From (2) £r (5) we obtain 


r dp 

’^'2az 

Now using (6) in equ. (4) 
we have - 

1(L^- 1= — 
pbaz dr 


dvp 

dr 


r i > 


\[ipj 


(^L- 

[dzl 


xp 


R<r<R, 


( 6 ) 


(7) 

( 8 ) 


dvc 


\x^\dz2 


,R<r<Rc 



-'^ = 0,0<r<Ry 
dr : 

Integrating equations (8), (9) and (10) 





Similarly 


1-/?^ +2ap(/? -1)| 


In non- dimensional form, we have 


s- + 2ad 


Where 


Volume flow rate is given by 


we 


Q = + Q 2 + Q 3 

R 

Q,=2n\v/dr^ 


( 17 ) 


Qj - 27t jV'/d/' and 


= 27 t \v/dr 

0 

Now, from equations (11), (12), (13) and (17) 
have 




3 r 


2ii, 


1/7" o 2 o h 2 
-- — + 2apa, ~2apha. 


2 2u2 2 U3^<^c 

-ac^ +ac^h^ +-^\ — --cL^h +-^ 
M'c [ 2 3 o 

Apparent -fluidityis given by 


F. = 


2Q 


-nt^Rp 


Using equation (18) in (19) 
we get 
_ 1 




■P L 


- - ^ + 2apac^ - 2apfiah - a/ + afh^ 
2 2 


( 18 ) 


(19) 


%cl2 3 




( 20 ) 
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Putting 






Resuit and Discussion - 

There are many different models of blood flow 
through narrow vessels with different boundary 
conditions and flow parameters but the scope of their 
quantitative study is very limited because of 
experimental problems. 

Here we have considered Up = 0 and ctp 0 i.e. 
peripheral layer is newtonian and non - Newtonian 
fluid respectively. 

Variation of apparent fluidity with peripheral layer 
thickness and yield stress (a^) are given below in tables 
1 and 2. Velocity idistributions are alos disscussed is 

Table ... 3 .. — • 

From fig 1 and 2, it is clear that apparent fluidity 
fa, increases with d (plasma layer thickness) and 

decreases with (yield stress). This indicates that 

when hematocrit decreases, apparent fluidity increases. 
The results are composed with results of Gupta £r 
Sehadri and found a good similarity in them.Also when 
'ap=0, fa is found. The values for Up =0 are greater 
than those for ap =0.01 and thus we;- found that 
apparent fluidity decreases with yield stress of plasma. 

Fig 3 shows that velocity field decreases fastly 
with r in the peripheral layer. In the core region. 



(Plasma layer thickness) 


Variation of Fa with 


Table 2 


(Varition of apparent fluidity with yield stress^ 


0.08 

0.10 

0.12 

0.14 

0.16 

0.18 

0.37 

0.421 

0.471 

0.522 

0.558 

0.624 

0.368 

0.420 

0.470 

0.520 

0.556 

0.623 




ai 6 . 

Fa . 

(a, = 0) 


0.559 ' 

Fa 

(a, = 0.1j 


0.558 













































Table - 3 

(h=40%, a/Rp=0.l, OL,=0.1, ap=0.01) 



(Variation of 


0 

0.2 

0.4 

0.6 

0.8 

1.0 

1 

0.96 

0.84 

0.64 

0.36 

0 

0.98 

0.94 

0.83 

0.63 

0.35 

0 





















Table - 4 


(Velocity Distribution in core region) 


r 

Rp 

0 

0.2 

0.4 

0.6 

0.8 

1.0 

h=.9 

Vc 

o(p=0.t) 

0.435 

0.45 7 

0.435 

0.3 70 

0.261 

0.103 

h =.9 

Vc 

c ^ p = 0.1 

0.433 

0.45 5 

0.43 3 

0.368 

0.259 

0.101 

/ 7=.85 

Vc 

a ( p = 0 

0.486 

0.507 

0.486 

0.421 

0.312 

0.159 

h =.85 

Vc 

<3?O=.01 

0.483 ■ 

0.504 

0.483 

0.418 

0.309 

0.146 


























veloity decreases slowly. On comparision of results with 
enperimental results of Bugliarello and sevilla we found 
a good similarity. When oc^ =0, velocity is greater then 
those obtained for ap=0.01. 
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interaction of blood and the vessels, hence one may 
extend the definition of apparent viscosity of steady 
flow to the oscillatory flow by replacing the steady 
Poiseuille flow by womersley's oscillatory flux and 
considering the equality of flow amplitudes. 

In present analysis we analyse the pulsatile blood 
flows throw straight and long narrow rigid circular 
tubes whose radius is constant, it is assumed that 
flow of blood is axially symmetric and in the z direction. 
The blood consist of two layers; one is a plasma of 
Newtonian fluid and the other is certral core of non 
Newtonian fluid. The thickness of Plasma layer is 
assumed to be constant and independent of time and 
location. 


6.2 The Governing Equaations - 

The momentum and continuity equations for fully 
developed viscous incomressible laminar flow under no 
body forces in the cylindrical coordinates Cr, 0, z), 
whose origin lies on the vessels axis are 


^ + l A| 

' dt ' dz r dr 
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where V is the axial component of velocity, P 
the peressure, t the time and p is the fluid density. 


The pressure gradient which is a function of time 
but not of axial loction, is, given by 


Where w is the angular frequency, p' and A are 


constant 


6.3 Constitutive Equations - 

Constitutive equation for the shear stress t and 
shear rate e is given by 


= Xg2 +}X^e^ 

where x/ is the yield stress, p the coefficient of 
viscosity. The equation is reduced to that for a 
Newtonian fluid when ^0=0 


The coeffi.cient of viscosi 
normal blood are given by 

u = uil + 0.025H + 7.35x10-'‘H^) 




Where M^p is the viscosity of plsma, H the 
hematocrit of blood and He is the hemotocrit 


concentration below which 


is negligible and it 


normally ranges from 4 to 8. The constant B is strong 
function of fibrinogen concentration and has value in 
the range 0.6 xlO’^ to 1.2x10’^ for normal blood. 

These equations are used in numerical results. 


6.4 The Boundary Conditions - 


s=ao-£ 






It is assumed that the fluid representing the blood 
■' contained .in a cylindripl tube of radius Ro- .consists 
of two 'layers - one representing the plasma near the 
tube wall having a constant thickness 5 and the other 
the central core. The inner radius of the tube has 
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Where and p are the densities and and a 
are the womersley frequency parametere of plasma and 
blooia respectively." Th^dimrence between V ahd p is 
so small that we put p/=p in the calculation. 

Using equations (12) and (13) in (11) we obtain 




been taken to be R = X Ro, where K = 


The boundary conditions are given by 


v, = o 

at /' = f?o 

(8) 

v,=v. 

at r = xRQ 

(9) 

T =^2 

at r = xRQ 

(10) 

II 

0 at r :=0 

. (1-1) 


The quantities in the plasma layer and in the 
central core are designated by the subscripts 1 and 2, 
respectively. 


6.5 The Mathematical Analysis - 

In the mathematical analysis we introduce the 
following non- dimensional variables and parameters. 


D . t = Wt 
^0 


( 12 ) 



the equations of velocity nd motion in plasma and 
core region. 


dVi 


dr 


dV, 


dr 


= ti 

^ J_ 

ho 




1\ 


T22 -(3? 

V y 


a' 


di 


-2P{t\ + --t{rT,] 
rdr^ > 


(14) 


(15) 


(16) 


By making use of the same equations (17) and 
(18), the boundary conditions are also transfromed 
into the following forms 


v; = o at Y=1 


at 


x^ = X 2 at r=x 


dV^ 

dr 


■ 0 at r = 


r = 0 


( 17 ) 


( 18 ) 


(19) 


-( 20 ) 


Now in the arterioles the womersley frequency 
parameter (a or a,) is very small compared with unity. 
Then the velocity \/ and shear stress x are assumed 
to possess expansions in tetris of the parameters a^ 
which are of the from 


V,=V^,+a2V, 


( 21 ) 


X/ = To/ +a^xi/ (2 2) 

Terms of order a”* of higher orders will be 

neglected. 


6.6 Velocity and shear stress - 

The solution of equations (14-16) with boundary 
conditions (17-19), gives the equations of velocity and 
shear stress distribution. 

Now using the equations (21) and (22) in (16) 
we get 


a 


dV oi 2 dV/' 
— =- + a — =- 

dt dt 




EquatitJ|. the coefficients of a° ^nd 


a 


2 we have 


0 = -2P|(| + 7|;f^») 


(24) 


dt r dr ^ ’ 


(25) 


integrating equs. (24) & (25) with boundary 

condition that shear stress is finite at the axis we 



r = T02 


(26) 


r uv 02 , ^ 

rxi 2 = r — ^dr + C. 


(28) 


Where and C 2 are constants 


Similarly, using equ's (21) and (22) in (14) and 
(15) and equating the coefficients of a° a^> we 


have 


(29) 


(30) 


dV 02 1 1 

8r 



equation C26) in equ. (29) and (31) and then intonating 
we have - 


we obtain 


From equations 


Where 


Now using equations (26), (35) and (36) in equatin 




X2=Pt]r 


1 + cx 


2j_^ 

P dt ‘ 


(3 8) 


where and L 2 are known functions of C} , 

}'), given by 
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7^0 


/ 1 l^ 

7x"-8y2x2 


(39) 


U=-{'\-x^] + z^\(2x^-A-~y^ 


1 


8|.io 


21 ' 


7x2 -4r2 


(40) 


Shear stress at the wall xw 's also an important 
quantity and is given by 


Xw 


= p(F) 


, 1 dp 


1 + a^-^4 

P dt , 


(41) 


Where 


4 


8 


(l-x"y+— 7x'-8y2x 
^ . f 7po V 


1 

/2 v2 


(42) 


Now for the velocity distribution, using equations 
(26), (35), ( 36 ) in (30), and applying boundary 

conditions ( 17 ), (18) we have- 


17. 

8 dt 


1-4-2 , - 3) 1 

-r -r *^Log,r*-]-— 


{ 17 


7x^-8y2x2 ^og^r 


(43) 



y 7x2r2-2r +5x 


From equations (21), (33), (34), (43) and (44) 


we obtain 


Where N^, N 2 and N 3 are know functions of 


given by 









6.7 Flow Rate - 

the volume flow rate Q (t) of the fluid flowing 
per unit time across the cross section of the vessel. 





is given by 


r^xRa 


r=Ro 


Q{t)= J 2%V2rdr+ jlnV/dr 


r=0 




4 Ol 
^l1 


7cRo"P 


1 

J V^rdr + jv^rdr 


( 50 ) 


Now, using equations (45) and (46) in (50) we 
obtain 



5 + a=J-®T 

! Qi vi 

p dt 


(51) 


Where 




Qi = — S and T are known functions of f , 


8|^i 


given by 


S = 1-X" + 


X' 


21n 


11 

21X-48y2x2 +28y 


0 ^ 


y 


(5 2) 


r=l(l-3X^ +2X®)+ ^ 


*) y 1 


7^0 V 


7x2 -8y2 


Log^X 


. . 479v^^^5 , A.,1 


+ • 


1^0 


6 462 35 J 


(53) 


The steady flow rate can be obained by substituting 
p (t) = l in equ. (51). 
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VPRIATrON OF reduced PRESSU/?E with REDUCti) 

flriPLxrupE 


VARrflTrON OF 1?&PUC&D PLUG] RflOXUS WCTW REDUCED 

BNpLrTUDE 




Results and Discussions 



The characteristies discussed are the pressure, 
shear sress at the wall, flow ^rate and visocsity. 

Here we have considred the case of vessel radius 
Rq = 10x10'^ cm, mean pressure gradient P’ = 0.62 x 10"^ 
dynes/cm^ of which plasma viscosity is p, =0.012 
centipose. 

The value of womrersely frequency parameter a 
is about 0.1 in arterioles. The experimentally measured 
values of viscosity of normal human blood varies from 
0.03 to 0.05 Poise. The value of yield stress for 
normal blood is of order 0.04 dynes per cm^. These 
values depend on various factors such as, hematocrit, 
fibrinogen concentration etc. 

It is seen that the effect of a on velocity, 
shear stress, plug radiug.-- and flow radius is very small 
in the range 0<a<1- 

. Fig. 6.1 shows the effect of reduced, amplitude A 
of pressure gradient on pulsatile blood flow. It is 

seen that’ the values of Rp, and Q0 v^ry 

greatlywith respect to time specially for large values 

of A. 




A 5a:/z 2 a 

VARtfiTroN OF R&DUC&J) HfiLL SH&AR STRESS V-lfTH 
RtOUClrD AHPLXTUOE 


vpRtArrON OF RtDuc&D flow Rate hetH. Reduced 

BMPirrooE 






The influence of reduced y\ield stress on R 


and Q\t) are shown in 


presence 


of plasma layer the value of xw not affected very 
much. The blood structure changes when the value of 
yield stress changes. The flow rate for steady flow 
can be seen from fig 6.2. It is seen that flow rate 
increases gradually with the increase in the thickness 

A 

of plasma layer. For smaller values of p the flow 


rate in central core is manimum and grdually decreases 
when the thickness of plasma layer is increased. 
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<lIhBOi'eticcil jVlodel of Steady ^lood J^low 
in /Harrow Vessel^ 



.4 


':j,y Ss^ts^odpiction : 

Mathematical modelling of blood flow in different vessels under-various 
physiological conditions have drawn the attention of several workers. Sin'ce 
poiseuilles empirical euqation of finding out the apparent viscosity of viscous 
fluid in fine glass capillaries, many workers attempted to find the apparent 
viscosity invivo and invitro, which is an important parameter affecting the blood 
rheology. Fahraeus and Lindquist found that apparent viscosity of blood 
decreases with the decrease of tub. diameter below 500 pm down to 7 pm. 
Recently, Kiani and Hudett (1991) developed the model assuming no yield stress 
behavio'ur of blood in these vessels. It has been recognised by several 
experimental observations and viscometric dates that blood possesses consid- 
enable amount of yield stress in narrow vessels. Scott Blair (1959), Charm and 
Kuriand (1968). bugliarello and Sevilla (1970), Oba (1981, 89) have proposeo 
casson constitutive equation to analyse the blood behaviour in narrow vessels 
at low and high she.ar rates. In the present analysis we have assumed that the 
fluid in core region satisfy casson equation and in marginal layer reg.on satis^ 
.ewtonian equation.' Apparent viscosity of blood has been determined as a^ 
function of yield stress, vessel diameter and peripheral layer (PPL) thickness. 














I 






rj 


1*2. J\flathamatical Analysis ; 

Incompressible fully developed laminar flow of blood is considered in a 
circular tube of small diameter. Two-layer flow model is assumed. Core fluid 

follow the casson equation. 




^ 1 


(1) 


and plasma peripheral layer follows the Newtonian equation 


P-p Y 


C2) 


Where -c is shear stress, Xy the yield stress, 7 the shear rates, casson 
core viscosity, jXp Newtonian plasma viscosity. 

Steady incompressible laminar flow equation - 

az r Sr 

after integration yields 


(3) 


'rz = 

2L 


(4) 


Where “ 

z, r - axial and radial measurements 

'^rz " shear stress normal to r in z direction 
p = pressure drop along length L of the tube. 

At the interfacial boundary between core and marginal fluids, shear 

stresses are equal and at the wall surface noslip condition are assumed. 
Substituting (1) and ( 2 ) in (4) we get 

CRw 2 - r2); R, < r ^ Rw 


Vp = P 


(5) 


4L p, 


V, 


-S- T, 


P , (Rc^ - rO - 4 


4Ljic 

(Rc - r) + 


3 

CRw^ 


f ? xyf 
2L/ 
Rc^) 


1 CR 3^2 . r3/2). 


C; 

Pc 

0 £ r s Rc 


( 6 ) 


4L^, 


Where Vp denotes plasma velocity and Vc denotes core velocity. 
Plug velocity is obtained by putting in equation (6) as 





+ P CRw^ - 


The total flow rate Q = Qp + Qc ■*■ *^1 

or Q = (\nrdr . V + r^‘2%rdr . V^ + it R,^ • V, 


r Q =-_^" M/if 

8L|I ' \ R.. l\P-r 


3 \ Rw'^ ^ ^ 

Where plug radius Rj = 'Idxy 

P 


The apparent blood viscosity M'app obtaine ^ ^ 

- /. . r/1 . '7RN4 .16 /I- ZSV^^ . 


+ 4 \ 

3 V. d / 21 

Where 


order smallness of 5. equation (9) reduces^ to 


Upto first 


Where FCa) = 1 ^ jf “ 

7 3 21 

X *1 o A» a 


,, . a ..asure^ent of /.e.d st«ss. th.fefor. at hlgV shear ta . a 

. chosen by Kianl and Hudetz 0991). For 

it tahes small value but not zero as chosen by it. 


constant F(a) and Ka), apparent viscosity decreases with the increase in 5, 
Fahraeus Lindquist effect can be observed from equation (10). Charm and 
ICurland have shown that 5 increases with decreasing a. We have taken the 
datas of Charm and ICurland and shown the variation of tigpp with a at constant 
d in table I. Assuming that the casson viscosity is related by general equation 
= Up (1-lCH) -1 

Where H is a function of hematocrit and K is an experimental constant, 
we can obtain at different concentration of hematocrit. 


apP 


<lA^£€ 9 

[Charm and Kurland (1968), Pp = .0125 poise, 
= .0259 "p^c. d = .01, H = .448] 


2.7790 

2.8293 

3.0584 


a 


F(a) 


f(a) 


.0005 

.0010 

.0050 


.9495 

.9291 

.8451 


.9559 
.9378 
, .8636 


2. 
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